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Abstract 

One of the main elements of the self-consistent generalized Langevin equation (SCGLE) theory 
of colloid dynamics [Phys. Rev. E 62, 3382 (2000); ibid 72, 031107 (2005)] is the introduction 
of exact short-time moment conditions in its formulation. The need to previously calculate these 
exact short-time properties constitutes a practical barrier for its application. In this note we report 
that a simplified version of this theory, in which this short-time information is eliminated, leads to 
the same results in the intermediate and long-time regimes. Deviations are only observed at short 
times, and are not qualitatively or quantitatively important. This is illustrated by comparing the 
two versions of the theory for representative model systems. 

PACS numbers: 64.70.Pf, 61.20.Gy, 47.57.J- 
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I. INTRODUCTION 



In recent work a new first-principles theory of dynamic arrest lias been proposed 



HQ. 



This consists essentially of the application o: 
tion (SCGLE) theory of colloid dynamics 



the self-consistent generalized Langevin equa- 



2, y, B, B 



to the description of the singular 



behavior characteristic of dynamic arrest phenomena in specific colloidal systems and con- 
ditions. The SCGLE theory was originally devised to describe tracer and collective diffusion 
properties of colloidal dispersions in the short- and intermediate-times regimes 0, B]. Its 
self-consistent character, however, introduces a non-linear dynamic feedback, leading to the 
prediction of dynamic arrest in these systems, similar to that exhibited by the mode cou- 
pling theory (MOT) of the ideal glass transition • The resulting theory of djTiamic arrest 
in colloidal dispersions was applied in recent work to describe the glass transition in three 
mono-disperse experimental model colloidal systems with specific (hard-sphere, screened 
electrostatic, and depletion) inter-particle effective forces The results indicate that 

the SCGLE theory of dynamic arrest has the same or better level of quantitative predictive 
power as conventional MCT, but is built on a completely independent conceptual basis, thus 
providing an alternative approach to the description of dynamic arrest phenomena. 

There is, however, a possible practical disadvantage of the SCGLE with respect to the 
MCT, and it refers to the fact that the MCT only requires the static structure factor of the 
system as an external input, whereas the SCGLE theory requires this information plus other 
additional static properties involved in the exact short-time conditions that the theory has 
built-in B. As it happens, however, the long-time asymptotic solutions of the relaxation 
equations that constitutes the SCGLE theory are independent of such exact short-time 
properties The questions then arise if a simplified version of the SCGLE theory, in 
which this short-time information is eliminated, could be proposed, and to what extent such 
a simpler theory will still provide a reliable representation of the dynamics of the colloidal 
system not only in the asymptotic long-time regime, but also at earlier stages. In what 
follows we demonstrate that there is a simple manner to build this simplified version of the 
SCGLE theory, and that it is virtually as accurate as the full version, even in the short- and 
intermediate-time regimes. This finding will greatly simplify the application of the SCGLE 
theory of dynamic arrest. 

Let us summarize the four distinct fundamental elements of the full self-consistent gen- 
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eralized Langevin equation theory of colloid dynamics. The first consists of general and 
exact memory-function expressions for the intermediate scattering function F{k, t) and its 
self component Fs{k,t), derived with the generalized Langevin equation (GLE) formalism 



10|, which in Laplace space read 

^ + 1+C{k,z) 

Fs{k,z) = ^— ^^ . (1.2) 

^ l+Cs{k,z) 

where Dq is the free- diffusion coefficient, S{k) is the static structure factor of the system, 
and C(fc, z) and Cs{k, z) are the corresponding memory functions. 

The second element is an approximate relationship between collective and self-dynamics. 
In the original proposal of the SCGLE theory [s], two possibilities, referred to as the ad- 
ditive and the multiplicative Vineyard-like approximations, were considered. The first ap- 
proximates the difference [C{k,t) — Cs(fc,t)], and the second the ratio [C{k,t) / Cs{k,t)\, of 
the memory functions, by their exact short-time limits, using the fact that the exact short- 
time expressions for these memory functions, denoted by C'^^^^jk, t) and C§^^^{k,t), are 
known in terms of equilibrium structural properties j^, The multiplicative approxima- 
tion was devised to describe more accurately the very early relaxation of F{k, t) [6], but the 



additive approximation was found to provide a more accurate prediction of dynamic arrest 
phenomena In this paper, for "full SCGLE theory" we refer to the theory that involves 
the additive Vineyard-hke approximation. 



C{k, t) = Cs{k, t) + t) - C|^^^(A;, t)]. (1.3) 

The third ingredient consists of the independent approximate determination of Fs{k,t) 
[or Cs{k,t)]. One intuitively expects that these /c-dependent self-diffusion properties should 
be simply related to the properties that describe the Brownian motion of individual particles, 
just like in the Gaussian approximation which expresses Fs{k,t) in terms of the mean- 



squared displacement (msd) (Ax(t))^ as Fs{k,t) = exp[—k'^{Ax(t)y/2]. We introduce an 
analogous approximate connection, but at the level of their respective memory functions. 



The memory function of (Ax(t)y is the so-called time- dependent friction function A({t). 
This function, normalized by the solvent friction (q, is the exact long wave-length limit of 
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Cs{k,t), i.e., limk^Q Cs{k,t) = A(*{t) = A((t)/(o- Thus, we interpolate Cs{k,t) between 
its two limits, namely, 

Cs{k,t) = Cl^''''{k,t) + [AC{t) - A(A;), (1.4) 

where 

x(k)^[i + {k/k,fr' (1.5) 



is a phenomenological interpolating 
that follows the main peak of S{k) 



unction, with kc being the position of the first minimum 
The fourth ingredient of our theory is another exact result, also derived within the GLE 

n 

approach p^, this time for A(*(t). This exact result may, upon a well-defined simplifying 
approximation, be converted into the following approximate but general expression [2i] 



AC(t) = — / dk 
^ ' 3(27r)'n ' 



F{k,t)Fs{k,t). (1.6) 



Sik) 

Eqs. f ll.ip -( fT!6l) constitute the full SCGLE theory of colloid dynamics. Besides the un- 
known dynamic properties, it involves the equilibrium properties S{k), C'^^^^{k,t) and 
C§^^^{k,t), determined by the methods of equilibrium statistical thermodynamics. We 
should also point out that Eqs. (11.11) and (11. 2p are exact results, and that Eq. (II. 6p derives 
from another exact result. Hence, it should not be a surprise that the same results are used 
by other theories; in fact, the same equations are employed in MCT. The difference lies, of 
course, in the the manner we relate and use them. In this sense, the distinctive elements of 
the SCGLE theory are the Vineyard-like approximation in Eq. (11.31) and the interpolating 
approximation in Eq. (11.41) . 

The simplified version of the SCGLE theory is now suggested by the form that these 
distinctive equations (Eqs. (11.31) and (II. 4p ) attain for times longer than the relaxation time 
of the functions C^^^^{k,t) and Cl^^^{k,t). Under those conditions, Eqs. f[L3|) and f[Li|l 
become, respectively, 

C{k,t) = Cs{k,t). (1.7) 

and 
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Cs{k,t) = [ACit)]m. (1.8) 

It is not difficult to see that the original self-consistent set of equations (involving Eqs. 
( 11. 3p and ( 11. 4p ) shares the same long-time asyrnptotic stationary solutions as its simplified 
version. Such stationary solutions are given by [2] 

and 

where 7 is the solution of the following equation 

7 6n^n Jo [X{k)S{k) + P7] [X{k) + F7] ' ^ ' ^ 

The parameter 7 is the long-time asymptotic value, of the msd, i.e., 7 = lim^^oo (Aa;(t))^. 
In the arrested states, this parameter is finite, representing the localization of the particles, 
whereas in the ergodic states it diverges. 

It is then natural to ask what the consequences would be of replacing Eqs. (11.31) and 
(ll.4p of the full SCGLE set of equations by the simpler approximations in Eqs. (11.71) and 
(II. 8p . that no longer contain the functions C''^^^^{k,t) and C§^^^{k,t). Our proposal of a 
simplified version of the SCGLE theory consists precisely of this replacement, so that the 
"simplified SCGLE theory" consists of the exact results in Eqs. (II. ip and ( 11. 2p along with 
Eqs. (II. 5p and (II. 6p . complemented by the closure approximations in Eqs. (I1.7p and (11.81) . 

We have made a systematic comparison of the various dynamic properties involved in the 
SCGLE theory, including the intermediate scattering function F{k,t), its self component 
Fs{k,t), and other tracer-diffusion properties such as the time-dependent friction function 
A(*{t), the mean squared displacement or the time-dependent diffusion coefficient D{t) = 



(Aa;(t))^/2t. As expected, the scenario of dynamic arrest exhibited by this simpler theory is 
identical to that provided by the full SCGLEscheme. This is probably not surprising since, 
as indicated above, both sets of dynamic equations share the same long-time asymptotic 
behavior and the same asymptotic stationary solutions. What is surprising, however, is 
the degree of accuracy of the simplified theory in the short- and intermediate-time regimes. 
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In our systematic comparison we considered systems with soft-sphere, hard-sphere, and 
repulsive Yukawa interactions, systems with attractive (Yukawa) interactions, systems in 
three and in two dimensions, and both, mono-disperse and bi-disperse systems, in all cases 
with similar conclusions, that we illustrate with the following examples. 

Thus, in Fig. [T] we plot the intermediate scattering function F{k, t) as a function of time, 
evaluated at the position kmax of the first maximum of the static structure factor of a soft- 
sphere system. The pair potential, in units of the thermal energy ksT = (3~^, is given by 
(3u{r) = ^^I'^Y" ~ {r%Y + 1 for < r < cr, and it vanishes for r > cr. The system in Fig. [1] 
corresponds to z/ = 18 and to the volume fractions = vmcr'^/G = 0.515, 0.612, and 0.613, 
and the static structure factor was calculated using the prescription of Verlet and Weis 12|. 
The heavy-solid lines correspond to the simplified version, and the heavy-dashed lines to the 
full version, of the SCGLE theory. The very first feature to notice is the virtual coincidence 
of the results of these two approximations; in fact, only for = 0.515 the difference is 
appreciable. For the other volume fractions the results are not distinguishable in the scale 
of the figure, and this includes the vicinity of the glass transition which, as can be seen in 
the figure, is predicted to occur at 0g = 0.613 for this soft-sphere system. For = 0.515 
we also show the Brownian dynamics data (solid circles) reported in Ref. [6], as well as the 
results of the SCGLE theory within the multiplicative Vineyard-like approximation (soft 
dashed lines), to recall the fact that the multiplicative approximation sometimes provides 
a slightly more accurate quantitative description of the initial relaxation of F{k, t) (see 
the inset). In a longer time-scale, as indicated in the main figure, the prediction of the 
overall relaxation provided by the SCGLE theory complemented with the additive and the 
multiplicative approximations is quite similar. Furthermore, as discussed in Ref. {2], the 
additive approximation provides a simpler and more accurate description of dynamic arrest, 
partly because these phenomena do not seem to depend strongly on the short-time behavior 
illustrated in the inset of Fig. [H Thus, from now on, we shall omit further reference to the 
multiplicative approximation. 

A similar situation is illustrated in Fig. [2], this time for a system of colloidal particles 
interacting though a hard-sphere potential of diameter a plus an additional long-ranged 
repulsive Yukawa tail of the form (3u{r) = K exp —z{r /a — l)/(r/cr), with z = 0.15 and 
K = 500. For a volume fraction 0i = 4.4 x 10~^, it corresponds to the conditions of Fig. 1 
of Ref. [1^. In our present figure, however, we compare the full and the simplified SCGLE 
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theory for the time- dependent diffusion coefficient D{t) and for the intermediate scattering 
functions F{k, t), F^{k, t), and F'^{k, t){= F{k, t) —F'^{k, t)) in the short-time regime, where 
these differences are expected to be larger. Our present comparison indicates that the new 
simphfied version of the SCGLE theory leads to essentially identical results, even in these 
time regimes. 

Just like the MOT has been extended to mixtures [3], the SCGLE has also been extended 
to multi-component colloidal systems Q, Also in this case the simplified version of the 
SCGLE theory provides virtually the same description as the full SCGLE scheme, but its 
practical application is far simpler. In Fig. [3] a comparison is presented for a binary Yukawa 
mixture, with pair potential [3uij{r) = ^jKiKj exp —z{r jo — \)/{r/a) (1 < i,j < 2) in 
which a fraction xi of the particles (species 1) interact with a charge parameter Ki = 100, 
the other fraction X2 (species 2) with K2 = 500. The volume fraction of the more interacting 
species is kept fixed at 02 = 2.2 x 10^^ and 0i takes the values 6.6 x 10~^ (right column), 
2.2 X 10~^ (middle column), and 7.25 x 10^^(left column), corresponding to xi = 0.75, 
xj = 0.5 and xi = 0.25. This figure corresponds to the same conditions as Fig. 3 of Ref. 



15l | and, as in our previous example, the simulated S{k) was employed. It simply confirms 
the general conclusions of this communication, namely, that the simplified SCGLE theory 
provides a description of the relaxation of concentration fluctuations in colloidal suspensions 
qualitatively and quantitatively virtually identical to the full SCGLE theory. Its practical 
implementation, however, is much simpler than either the full SCGLE or the MCT schemes. 
This has simplified th e ap plication of the SCGLE theory to the discussion of dynamic arrest 



in colloidal mixtures [16| and in colloidal fluids adsorbed in model porous media [17|, that 



we report in separate communications. 
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FIGURES AND CAPTIONS: 




FIG. 1: F{k, t) as a function of time (in units of to = o'^Z-^o) at the position kmax of the first max- 
imum of S{k) of a soft-sphere system with z/ = 18 and volume fractions (j) = 0.515, 0.56, 0.60, 0.612, 
and 0.613, calculated with the simplified (heavy solid line) and the full (heavy dashed line) versions 
of the SCGLE theory. For (p = 0.515 we also show the intermediate-time Brownian dynamics data 
(solid circles) and the results of the SCGLE theory with the multiplicative approximation reported 



in Ref. 



The inset is a close up view in linear scale of the short-time relaxation of F{k,t). 
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FIG. 2: Time-dependent diffusion coefficient D(t) and (self, distinct, and total) intermediate scat- 
tering functions F'^{k, t), F'^{k, t) and F{k, i) of a repulsive Yukawa system with z = 0.15, K = 500, 
and (f) = 4.4 x 10"''. Results of the simplified (solid lines) and of the full (dashed lines) SCGLE the- 
ory for the intermediate scattering functions are shown for times t = to (upper or more structured 
curves) and t = Wto. The solid circles are the Brownian dynamic data of Ref. [l^ . 



(|),=7.25x10'^ <f =2.2x10''^ |^=6.6x^6 




FIG. 3: Time-dependent diffusion coefficients Di{t) and D2{t) as a function of time (upper curve 
corresponding to species 1), and total intermediate scattering functions Fii{k,t) and F22{k,t) of 
a repulsive Yukawa mixture with z = 0.15, Ki = 100, K2 = 500 for t = tQ and t = lOto- The 
volume fraction of the more interacting species is kept fixed at (/>2 = 2.2 x 10~^ and cpi takes the 
values (f)i = 7.25 x 10~^ (left column), 2.2 x 10~^(center column), and 6.6 x 10~^(rigth column). 
The symbology and conventions are the same as in Fig. [2j 
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